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Abstract. In order to characterise the behaviour of 
Water Cherenkov Detectors (WCD) under a sudden 
increase of 1 GeV - 1 TeV background photons from 
a Gamma Ray Burst (GRB), simulations were con- 
ducted and compared to data acquired by the WCD 
of the Large Aperture GRB Observatory (LAGO). 
The LAGO operates arrays of WCD at high altitude 
to detect GRBs using the single particle technique. 
The LAGO sensitivity to GRBs is derived from the 
reported simulations of the gamma initiated particle 
showers in the atmosphere and the WCD response 
to secondaries. 
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I. Introduction 

Gamma Ray Bursts (GRBs) have captured the atten- 
tion of the asttophysics community during the second 
half of the past century, and are nowadays still a hot 
topic of study. Instruments in satellites and ground based 
experiments around the planet have been designed in 
order to study this impressive phenomena (see [1], [2] 
and references therein). Large data set of GRB properties 
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can be found in [2] for example. 

GRB emissions in the keV to MeV energy range has 
been studied in great detail. Though high energy GRB 
emissions have been reported (for example up to 2 
GeV for GRB940217 seen by EGRET), there is little 
understanding on the subject. Satellite Fermi will be in 
charge of this mission and ground based experiments 
will also conttibute. Theoretical predictions on the high 
energy emission in relativistic fireball models are given 
in [3], [4], [5]. Feasibility of ground based observation 
of Gamma Ray Bursts in the lGeV and 1 TeV energy 
range using the singe particle technique has been exten- 
sively analysed (see for example [6], [7], [8], [9]). The 
Large aperture Gamma-ray-bursts Observatory (LAGO) 
is designed to observe high energy GRB emissions by 
the single particle technique using water Cherenkov 
detectors (WCD) (see [10] in this volume). It consists 
of various sites, located around the world, with large 
GRB detection efficiency, either for their size or their 
altitude. The LAGO has 2 sites in operation, they are: 
Siem Negra (4557 m a.s.l.), Mexico and Chacaltaya 
(5200 m a.s.l.), Bolivia. One site at Pico Espejo (4765 
m a.s.l.), Venezuela is at present under construction and 
new sites are now under consideration in Argentina, 
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Peru, Colombia, Guatemala and Himalaya. High energy 
photons from a GRB are expected to produce particle 
cascades after reaching the atmosphere, thought with few 
secondary particles reaching the Earth surface. Neverthe- 
less, many photons detectable by WCD are expected to 
arrive during a GRB. 

In this paper we study the detection prospect of a GRB at 
three high altitude LAGO sites: Sierra Negra, Chacaltaya 
and Pico Espejo. In order to achieve this end, we 
simulated billions of showers initiated by photons with 
energies between 0.5 GeV and 1 TeV generated above 
each site, and the response of WCDs to secondaries 
reaching the ground level. 

II. Photon Showers at Ground Level 

In this section we investigate photon shower at ground 
level for three sites in the LAGO Observatory, Cha- 
caltaya (16°20' S, 68° 10' W , 5200 m a.s.l.), Bolivia; 
Sierra Negra (18°59' N, 97° 18' W , 4557 m a.s.l.), 
Mexico and Pico Espejo (8°38' N, 71°09' W , 4765 m 
a.s.l.), Venezuela. Primary energies are considered from 
0.5 GeV to 1 TeV with zenith angles between 0° to 
30°. Showers were simulated using the Extensive Air 
Showers package AIRES[11]. No thinnig was applied, 
thus all particles were followed until absorbed in the 
atmosphere or reached the ground. We simulated show- 
ers with and without geomagnetic field and found no 
significant difference in the shower development at these 
latitudes. We analyse the nature, number and energy 
spectrum of secondary particles reaching the ground 
level as function of the photon primary energy and zenith 
angle of injection. 
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In figure la we represent the mean number of particles 
(photons, e + /e~ and /i + //i~) reaching the ground level 
as function of the primary photon energy at the three 
sites. The number of particles at ground level increases 
with the primary energy as a power law for the sites 
with index shown in table 1. The mean number of 
secondary photons are about one order of magnitude 
above the mean number of e + /e _ , while the number 
of muons start to be important after ~ 1 TeV for 
the three sites. Already some particles arrive for the 
lowest primary energies while hundreds and thousands 
of secondary photons arrive around 100 GeV and 300 
GeV respectively (see table II). This is an advantage for 
surface detectors able to catch secondary photons. 
We compare the development of vertical showers with 
the ones slanted at 30° of zenith angles. The mean 
number of all secondary particles reaching the ground 
level vs primary energy for these two angles of injection 



at the three sites are displayed in figure lb. Significant 
differences are observed for the three sites, which means 
that high altitude surface detectors are very sensible to 
the primary particle injection angle with respect to the 
detectors in the atmosphere. Next section is dedicated to 
the response of the surface detectors to vertical showers, 
studies depending on the zenith angle will be reported 
else where[12] 

III. Water Cherenkov detectors response to 

HIGH ENERGY PHOTON SHOWERS 

The LAGO aims to detect secondary particles from 
showers generated by photons using the single particle 
technique [7] which essentially involves the detection 
of secondary particle by individual WCDs. As photon 
showers around GeV give only isolated particles at 
ground level, even at high altitudes, it is not easy to 
reconstruct the shower development. Nevertheless, if 
secondary particles generated by these gamma rays rep- 
resent a significant excess of events on the background 
a GRB would be detected in coincidence with satellites. 
The advantage of WCD relies in the fact that during a 
GRB many photons are expect to arrive, most of them 
able to create pairs in the water tanks with sufficient 
energy for Cherenkov emission. 

From our previous study on the features of the particles 
at ground level at each site, we are now able to analyze 
the response of a WCD to these secondary particles. 
Each site consists of 3 x(~ 4 m 2 ) WCDs placed in 
an irregular array (there are no specifications in the 
separation of the tanks since we are using the single 
particle technique). For the case of these three LAGO 
sites, the efficiency to GRB detection is due to the high 
altitude of the sites [8] (details about the experimental 
set up can be found in [13]). 

We simulated the WCD response to secondary particles 
reaching the ground level, produced in showers initiated 
by high energy photons, for each high altitude LAGO 
site. In order to achieve this end, we used a fast WCDs 
simulator developed by the LAGO collaboration. The 
main parameters considered for the detectors are the 
tanks geometry, the water altitude; and the number, 
geometry and position of PMTs. We simulated photons, 
e + /e _ and entering vertically in the tanks 1 , for 

energies between 10 MeV and 10 GeV : . Energy spec- 
trum of secondary particles reaching the ground level for 
the three sites, produced in electromagnetic showers with 
10 GeV and 100 GeV of primary energy is presented in 
figure 2a as an example. Similar behavior is found for the 
three sites, while Chacaltaya highlights over the others 
two. For both primary energies we observe that most of 
the particles arrive with energies about 10 MeV, then 
for E — 10 GeV the spectrum presents a steepening 

'It has been reported in [9] that the response of WCD does not 
depend much on the incidence angle of the secondary particles except 
at very low energies such that electrons can be absorbed by the plastic 
layer of the detector. 

2 AIRES output present histograms from 10 MeV to 7079 GeV 
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Fig. 1. Left: Mean number of secondary photons (black), e + /e _ (white) and /i + /fi~ (gray), reaching the groung level vs primary energy. 
Right: Mean number of all secondary particles reaching the ground level vs primary energy for zenith angles 0° (black) and 30° (white). In both 
cases circles, triangles and squares stand for Chacaltaya, Pico Espejo and Sierra Negra respectively. 
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Fig. 2. Left: Energy spectrum of secondary particles reaching the ground level al Chacaltaya(circles), Pico Espejo (triangles) and Sierra Negra 
(squares) for primary energies 10 GeV (white) and 100 GeV (black). Right: Probability of particle detection at a 4m 2 WCD for energy range 
from 10 MeV to 40 GeV 



after ~ 14 MeV and almost fade away over 0.2 GeV 
whereas for E — 100 GeV the steepening appears after 
~ 20 MeV and we can still see particles over the 1 GeV. 
The threshold energy for e + /e _ Cherenkov production 
in pure water is ~ 0.77 MeV making photons with 
E 1 > 1.55 MeV suitable to crate pairs detectable for 
our PMTs, thus our WCD response simulations are good 
enough in the energy spectrum of secondaries obtained 
with shower simulations. 

Detection probability P(E S ) for each 4m 2 tank is plot- 
ted in figure 2b which is defined as 

where n tr i g (E s ) is the total number of triggers and 
N tota i(E s ) is the total number of particles entering in 
the tank. We used 4 ADC as threshold for our trigger 
criteria. In this plot we can appreciate that detection 
probability rapidly increase between E s = 10 MeV and 
E s = 100 and saturates near ~ 1000 MeV 
We calculated the mean number of triggers as function 



of the primary energy from the convolution of detection 
probability of particles hitting the WCD, and the differ- 
ential spectrum of secondary particles, for each primary 
energy, taken from the shower simulations. Results for 
vertical showers are shown in figure 3a. In this plot we 
can observe an increment of the mean number of triggers 
with the primary energy, being Chacaltaya clearly more 
sensitive to low energy photons than Pico Espejo and 
Sierra Negra (both overlapped along their curves). 

IV. GRB DETECTION AT HIGH ALTITUDE LAGO 
SITES 

The LAGO is designed to detect GRBs by a sudden 
increase of the background in a short period of time. 
That means the counting rate in a surface observatory, 
summed over all the tanks, should be significantly higher 
than the background fluctuations. In this section we 
study the probability of GRB detection by arrays of 
WCD placed at Chacaltaya, Pico Espejo and Sierra 
Negra individually. The flux of a GRB at Earth can be 
described by the power law (See 2006 BATSE catalog 
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Fig. 3. Left: Mean number of triggers as function of the primary energy for Chacaltaya (circles), Pico Espejo (triangles) and Sierra Negra 
(squares). Right: Minimal fluences vs the spectral parameter a for different E max at each site, black symbols stand for Chacaltaya, white for 
Pico Espejo and Gray for Sierra Neg Fa. J^max — 

lGeV circles, E max = 10 GeV squares, Ernax — 100 GeV up triangles and Emax — 1 

TeV down triangles 



[2]) 



J{E) = K 



(I) 



(2) 



where if is a spectrum normalization constant, Eq is 
fixed to lGeV and the spectral index a takes values 
between 2 and 3. The mean value of a is around 2.2 
and no high energy cut-off has been observed [14], [15]. 
We will naively assume that such power law is still 
valid for high energy photons. We compute the mean 
number of triggered particles generated during a GRB 
integrating the typical GRB flux and the mean number of 
triggered particles from a shower generated by a photon 
with energy (< N(E) tr i g >)over the primary energies. 

,iTeV 

< ntrigger(a) >= / J(E) < n tna {E) > dE . (3) 

Jo.sGeV 

In this calculation, K is still a free parameter. Minimal 
K can be calculated from the condition: 



< n tr igg e r(a) > > 



B 



(9) A V AtN d ' 



(4) 



here s is the statistical significance, 9 is the zenith angle, 
B is the background counting rate, At is the duration 
of the flux of photons, A is the area of the WCD and 
Nd is the number of WCDs at each site. We will take 
At — Is, s = 5, as we are working with vertical showers 
for the tanks response then 9 = 0, B = 12 kHz per 
tank (which is the background at Chacaltaya), our tanks 
have typically A = A m 2 and Nd = 3. Once minimal 
K(a) is calculated, sweeping on various values of the 
spectral parameter, then one can estimate the properties 
of a possible GRB which can be detected by our surface 
detectors. The GRB flux at Earth can be calculated for 
different values of a from its luminosity L at certain 
redshift z in a range of energies and the distance from 
the source r(z) as follows 



J(E)EdE . 



(5) 



An estimation about the minimal fluences that we could 
expect to detect with our surface detectors can be calcu- 
lated by the integral in 5 from 0.5 GeV to E max , where 
E m ax varies from 1 GeV to 1 TeV. Results of these GRB 
fluences as function of the spectral parameter a are rep- 
resented in figure 3b for E max = 1, 10 , 100 and 1000 



GeV. We used At = 1 s. We take z<l supposing small 
distortions in the spectrum. Fluences shown in figure 3b 
exhibit that our surface detectors seems to be sensitive 
in the range of 10~ 7 and 10~ 6 erg/cm 2 . 
These results are comparable with the prediction made 
in [7] for high altitude (4.5 Km a.s.l. to 5.5 Km a.s.l.) 
surface detectors of 10 3 m 2 scintillators array. It is 
reasonable to expect sensitivity to even smaller fluences 
since we use WCDs able to detect the huge amount of 
secondary photons coming from high energy primary 
photons. The Pierre Auger Observatory shows sensitivity 
to high energy (between 10 MeV to 1 TeV) GRB photon 
fluences, using the single particle technique, between 
10~ 6 and 10" 3 erg/cm 2 [16]. BATSE detectors were 



reported[2] to be sensitive to fluences of 



io- 



the rage of energies from 25KeV to 2MeV among 2704 
events identified as GRBs. 

V. Discussion and conclusions 

Our atmospheric cascade simulations showed signifi- 
cant differences for zenith angles between 0° and 30°. 
Since high altitude surface detectors would be very 
sensitive to injection angle, inclusion of slanted showers 
in the response of the WCD and the extension of this 
study to larger zenith angles (> 30°) would be relevant. 
On the other hand, shower simulation results were not 
sensible to the geomagnetic field for the sites considered 
(with different latitudes). Nevertheless, this has to be 
checked for greater zenith angles [12]. 
Detection of GRBs with WCDs is powered by the 
presence of high energy pairs e + je~ into the water, 
created by photons (E 1 > 1.55 MeV), able to pro- 
duce Cherenkov light. Our results showed a significant 
amount of particles arriving with energies close to 10 
MeV even for primary energies of 1 GeV. Since the 
number of photons arriving to the detectors is about 
one order of magnitude higher than charged particles, 
their contribution to the signal is valuable. Moreover, 
the three sites studied in this work, Chacaltaya, Sierra 
Negra and Pico Espejo have a notable potential on 
GRB high energy photons detection, using the single 
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particle technique, due its high altitude. In the group, 
Chacaltaya (5200 m a.s.l., ~ 640 m higher than Sierra 
Negra and ~ 440 m higher than Pico Espejo) highlights 
with a minimal fluence lower than Pico Espejo and 
Sierra Negra by an order of magnitude, meanwhile 
detectable fluences appear very close for the last two 
(for example, Pico Espejo is 0.1 lower than Sierra Negra 
for E max = 1 GeV and a = 2.3). Our analisys showed 
that the LAGO high altitude sites (with 3 tanks at each 
site) are as or more sensitive for GRB detection than the 
complete Pierre Auger Observatory (1600 tanks) 
This results will allow us to analyze the space of param- 
eters involved in GRB detection, using the single particle 
technique, and propose limits in order to compare with 
satellite observations [10] 
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